How to realize a robust practical Majorana chain in a quantum dot-superconductor 

linear array 
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Semiconducting nanowires in proximity to superconductors are promising experimental systems 
for Majorana fermions, which may uhimately be used as building blocks for topological quantum 
computers. A serious challenge in the experimental realization of the Majorana fermions is the 
supression of topological superconductivity by disorder. We show that Majorana fermions pro- 
tected by a robust topological gap can occur at the ends of a chain of quantum dots connected 
by s-wave superconductors. In the appropriate parameter regime, we establish that the quantum 
dot/superconductor system is equivalent to a ID Kitaev chain, which can be tuned to be in a ro- 
bust topological phase with Majorana end modes even in the case where the quantum dots and 
superconductors are both strongly disordered. Such a spin-orbit coupled quantum dot - s-wave 
superconductor array provides an ideal experimental platform for the observation of non-Abelian 
Majorana modes. 

PACS numbers: 



Introduction: Solid state Majorana fermions (MFs) 
with non-Abelian statistics are the most promising can- 
didates for realizing systems with a topologically degen- 
erate ground state. Such particles have attracted recent 
attention both due to their fundamental interest as a new 
type of particle with non-Abelian statistics and their po- 
tential application in topological quantum computation 
[ll-Q]. Over the past few years it has been realized that 
topological superconductors provide one of the concep- 
tually simplest platforms for the solid state realization of 
MFs. More recently it was shown that topological super- 
conductors with MFs can occur fairly generically in sys- 
tems containing the three ingredients of conventional s- 
wave superconductivity, time-reversal breaking (e.g. via 
a magnetic field), and spin-orbit coupling [5]. A sim- 
ple topological superconducting (TS) system supporting 
MFs, which has attracted serious experimental attention 
consists of a semiconductor nanowire in a magnetic 
field placed on an ordinary superconductor Q. The s- 
wave proximity effect on a InAs quantum wire, which 
also has a sizable SO coupling, has already been realized 
in experiments [7j, and the search for the corresponding 
TS phase, in the presence of a magnetic field, with local- 
ized Majorana modes in semiconductors is actively being 
pursued in many laboratories worldwide. The theoretical 
demonstration 9] that such ID systems can show non- 
Abelian statistics has made these systems viable compo- 
nents for quantum information processing. 

While the proposal for the realization of MFs in semi- 
conductor structures has generated considerable theoret- 
ical interest and even inspired significant experimental 
efforts Q, a few key challenges relating to non-idealities 
in realistic systems remain. One such practical challenge 
is whether it is possible to control the chemical potential 
of a semiconductor wire, which is absolutely essential to 



tune the system into the TS phase, since it is in contact 
with a superconductor. The other key challenge arises 
from the fact that the TS gap is typically suppressed by 
strong disorder .1C| . Therefore the present proposals for 
realizing MFs would either require clean semiconductor 
systems or a spin-orbit coupling which is large enough 
to dominate over the disorder scattering effects ^lOi]. It 
is completely unclear at this stage whether these key 
obstacles of chemical potential tuning and low disorder 
can be overcome in the currently studied semiconductor 
nanowire structures unless substantial materials develop- 
ment occurs first. Our proposal in this paper overcomes 
both of these problems completely in one stroke by com- 
ing up with a new architecture (i.e. topology) for creating 
the Majorana- additionally, our proposal involves using 
semiconductor quantum dots coupled to superconduct- 
ing grains which is in many ways a simpler system to 
work with than the earlier two-dimensional 5] or one- 
dimensional structures Q for observing the TS phase. 
We are fairly confident that the system proposed in the 
current paper is by far the most ideal experimental plat- 
form to search for the solid state Majorana modes. 

In this paper we describe a new system composed of 
a linear quantum-dot(QD)/superconductor (SC) array 
(shown in Fig. [Ua)) to realize non-Abelian Majorana 
fermion modes at the ends of the array. The effects of 
disorder in the QD/SC array are suppressed in a trans- 
parent way because the topological state can be formed 
as long as the electron can hop from one QD to a neigh- 
boring QD in the array, which is a much easier condi- 
tion to satisfy than the near ballistic transport condition 
necessary through a wire. Strong disorder can typically 
renormalize the hopping across QD or SC islands, how- 
ever as long as the disorder energy scale is less than the 
fermi energy, one can choose the QD/SC to be shorter 
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FIG. 1: (a) Geometry to create disorder robust TS systems 
with end MFs using a QD-SC array. In the appropriate pa- 
rameter regime non-Abelian MFs are predicted to occur and 
be localized at the end QDs of a disordered QD/SC array. 
Each QD/SC island in the array has a different index to indi- 
cate that each QD or SC is can have a different local disorder 
potential, hopping amplitude or proximity-induced supercon- 
ducting pairing potential. The chemical potential fi„ in each 
QD is controlled by the gate„. Tunneling through the super- 
conductor is assumed to be weaker than the magnetic field 
induced Zeeman splitting SVz,n ~ Vz of the QD levels. Ap- 
plying fluxes to the SC islands can compensate relative signs 
between anomalous and normal tunneling, (b) Level structure 
of the n*^ QD relative to the fermi level (shown by a dashed 
line). The level spacing 5 En is assumed to be the largest en- 
ergy scale (~ 2 — 3 meV), the spin splitting SVz,n is of the 
order ~ 0.3 meV. The spin texture in each level, which is a 
result of spin-orbit coupling and is determined by the length 
of the QD relative to the spin-orbit length Iso (see Eq. [3|, 
plays a crucial role in allowing spin-singlet proximity-effect 
between different QD„ and QDn+i. This level structure is 
qualitatively independent of local disorder potential. 



than the locahzation length to enable significant trans- 
mission probability. The confinement potential of the 
QD separates the electronic energy levels in the dots as 
shown in Fig. [Hb), so that there are two levels (because 
of time-reversal symmetry induced Kramers-degeneracy) . 
For clean TS systems, it is known that gapped super- 
conductors are topological as long as they have an odd 
number of occupied bands at the fermi-level [1, 11|. To 
realize a TS QD/SC array, the Kramer-degeneracy of the 
levels must be broken so that an odd number of occupied 
levels per QD occurs at the fermi level. The Kramer's de- 
generacy can be lifted as shown in Fig. [IJb) by breaking 
time-reversal symmetry by an externally applied mag- 
netic field-induced Zeeman splitting Vz- Disorder in the 
QD affects the wave-function and causes fluctuations in 
the spacing of the levels in the dots, but does not af- 
fect the qualitative structure shown in Fig. [IJb). In the 
limit where the QD/SC tunneling t (which is physically 



characterized by the inverse tunneling time from QD to 
SC) is small compared to the Zeeman splitting between 
the spin-split levels, the chemical potential can be ad- 
justed by using gate voltages (shown in Fig. [Ha)) so that 
only a single level in each QD participates in transport 
through the QD/SC array structure. As in the TS phase 
of other [HI, Q superconductor-semiconductor structures, 
spin-orbit coupling plays a crucial role in determining 
the superconducting gap in the TS phase. Because the 
single levels in each QD are spin-polarized almost in the 
same direction, spin-orbit coupling is necessary to lead 
to a rotation of the spin-polarization between the left 
and right end of the QD (as shown in Fig. [Ijb)) so as to 
allow a proximity-induced SC coupling between neigh- 
boring QDs. Note that the on-site proximity-induced SC 
coupling in the same QD is forbidden by the Pauli exclu- 
sion principle. Furthermore, Coulomb charging effects 
can be ignored in the QD since we are operating in the 
parameter regime containing a single active level per QD. 
We mention that the number of electrons per dot (they 
could be as few as a few 10s), whether the dots are iden- 
tical or not, whether hopping through the whole linear 
array is coherent or not, how many dots there are in the 
system (as few as 5-10 should work) are all non-issues 
for the realization of the TS phase (as we have explic- 
itly verified numerically)- the only constraints are that 
there should be only one electron at the chemical poten- 
tial in each dot (which can be secured by adjusting the 
local gates individually on each dot) and there should be 
enough hopping to ensure local proximity effect. 

While for perfectly periodic QD/SC arrays, the condi- 
tion for realizing a TS systems is given by whether there 
are an odd number of partially filled levels in each QD, 
one must understand the Bogoliubov-de Gennes (BdG) 
Hamiltonian of the QD/SC array more generally to un- 
derstand the condition for the emergence of MFs in the 
disordered QD/SC array. In the limit of weak (compared 
to level spacing in the QD) QD/SC tunneling, the low- 
energy part of the BdG quasiparticle spectrum in this 
system, containing the zero-energy Majorana modes in 
the TS phase, can be described in terms of an effective 
ID Kitaev chain ll| lattice-model 
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where ipn is the creation operator in the active level in 
each QD, fin the effective chemical potential in each QD, 
tn and A„ are the normal and anomalous tunneling am- 
plitudes between the QDs. Note that in a disordered 
QD/SC array, the amplitudes t„,A„ would vary along 
the length of the chain. The complex conjugation sym- 
metry of the full BdG Hamiltonian for the QD/SC system 
will be shown to typically result in real effective param- 
eters i„, A„. 

The effective Kitaev chain description of the QD/SC 



3 



array allows us to reduce the complexity of an array 
of disordered QDs and SCs to just three parameters 
^„,t„, A„ for each QD. The energy of the QD level 
relative to the Fermi level can be controlled by a gate 
voltage so that only a single level in each QD is par- 
tially occupied. In the perfectly periodic case where the 
amplitudes i„ and A„ and the local potential /i„ are 
independent of n, the topological condition is given by 
l^ril < Knlill]- In a disordered QD/SC array, the magni- 
tudes of tn can be controlled to a certain degree by con- 
trolling the QD/SC tunneling. However, the magnitude 
of A„ or the signs of tn, A„ are more difficult to control 
given the disorder in the SC islands and the short fermi 
wave- vectors of electrons in the SC islands. In principle, 
such a strongly disordered Kitaev chain is expected to 
have a large number of low-energy states, similar to a dis- 
ordered topological superconductor, hindering the obser- 
vation of non-Abelian statistics and consequently topo- 
logical quantum computation. We will show that these 
low-energy states can be eliminated by shifting the phase 
of the superconducting order parameters on a sub-set of 
the SC islands by tt. By applying the appropriate phase 
shifts a robust TS phase with end MFs and a bulk gap 
of order ~ A„ with MFs at the ends can be obtained in 
the disordered ID QD/SC array whenever the conditions 
and sign(t„A„) = sign(i„+i A„+i) are satisfied. 
More generally, we show in the appendix that localized 
MFs are obtained whenever the localization length of an 
effective Hamiltonian is longer than the superconducting 
coherence length. Since A„ and i„ are of the same order 
this turns out to be easy to satisfy as demonstrated by 
the exact diagonalization of large number of short Kitaev 
chains. Phase shifts applied to the SC islands by values 
other than tt allow us to transform the Kitaev chain to 
this form even for complex hopping tn and pairing A„. 
Any segment of the linear QD/SC array can be tuned 
to the non-topological phase by tuning the chemical po- 
tential in the segment so that |At„+i| > max(|t„|, |A„|). 
Therefore we will establish that a gate-tunable robust 
TS phase can be obtained in a disordered linear QD/SC 
array as long as disorder does not suppress any of the 
tunneling parameters tn to being vanishingly small. The 
TS phase would have a robust gap that would be deter- 
mined by A„ with non-Abelian MFs at the ends of the 
QD/SC array, thus making this linear QD/SC array an 
ideal practical platform for the realization of non-Abelian 
statistics and MFs. 



Spin-orbit coupled quantum dots To derive the effec- 
tive Kitaev chain model given in Eq. [1] for the QD/SC 
array one must first understand the interplay of confine- 
ment, Zeeman splitting, spin-orbit coupling, and disorder 
in determining the wave-function of the active level in 
the QD in the vicinity of the fermi-level. Such an under- 
standing is obtained by considering a Rashba spin-orbit 



coupled ID system with a Hamiltonian 

[-^dl + laayd^ + +Vz(J, + F(a;)]^(:E) = E^jix) (2) 
m* 

where m* ^ 0.04me is the effective mass in InAs, Vz 
is the magnetic-field induced Zeeman potential, V{x) 
represents the combination of the confinement potential 
and local disorder potential and a is the strength of the 
Rashba spin-orbit coupling. We have approximated the 
QD to be one-dimensional since the confiment energy of 
InAs QDs along the transverse direction can be as large 
as ^ 40 meV [l^- The Rashba spin-orbit coupling leads 
to a variation of the local spin-density along the axis x 
of the QD, which can be described by transforming each 
eigenstate according to a spin-dependent gauge transfor- 
mation 

i;{x) '4>{x) = e''"'y''^'^°iP{x), (3) 

where Iso = :;^r^ is the effective spin-orbit length, which 
is the length travelled by an electron before it precesses 
by TT under the influence of the spin-orbit magnetic field. 
Substituting this into the Schrodinger equation we find 
that the spin-orbit term can be completely cancelled in 
this new basis so that the new Schrodinger equation is 
written as 

[ — + Vz{<7. cos + CT^ sin + V{x)]^{x) 
m* Iso ISO 

= iE + -^mx). (4) 

Ignoring the relatively small Zeeman splitting Vz ~ 0.5 
meV, the energy splitting for levels in the QD 5E ~ 
where vp ^ ^j2E/m* is the mean fermi veloctiy at an 
energy E relative to the bottom of the band and 2a is 
the length of the QD. Since the energy E must be less 
than the inter sub-band separation ~ 40 meV, we choose 
E ^ 5 meV, so that SE ^ 3 meV for a QD of length 2a ~ 
120 nm. In the absence of a Zeeman splitting (i.e. Vz = 
0), the energy levels in the QD, which are separated by 
around SE ~ 3 meV on average, are two-fold degenerate 
due to time-reversal symmetry. Disorder potentials can 
lead to statistical variations in the spacings of the QD 
levels. However, by tuning the gate voltage of the QD 
it should be possible to tune to a pair of levels that are 
separated from other levels by the average separation of 
6E ~ 3 meV. While the disorder potential in the QD 
is allowed to be large, we will assume that the length of 
the QD, 2a, is not significantly longer than the mean-free 
path so that the spatial part of the wave-function ^'0(2^) 
of the pair of levels in the QD extends over the entire 
length of the QD. 

Since Vz ^ 0.5 meV <^ 6E ^ 3 meV, the effect of 
Vz can be treated as a perturbation splitting the two- 
fold degeneracy of the pair of states at the fermi level. 
For a < IsOy the Zeeman field perturbation propor- 
tional to Vz in Eq. U in the degenerate space is given 
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by an eflFective magnetic field Vz[(Jz {ijjo\cosj^ Itjjo) + 
o'x {ipo\ sin \^po)] for the pair of states at the fermi 
level. The Zeeman field splits the degenerate pair of 
levels into non-degenerate levels separated by SVz = 

Vz^\ (^ol cos l^o) P + I (^ol sin ^ jV'o) P, which we 

assume to be larger than 6Vz > ~ 0.25 meV for 
a < I so- Moreover the wave- function of the lower of the 
spin-split states is ip{x) — ^po{x)e'^^'^y \az = —1) where 
tan 77 = (V'ol cos \ipo) / (V-ol sin j^\ipo)- Substituting 
back the wave-function ^'(2;) into Eq.|3l we obtain the ap- 
proximate expression for the wave-function of the active 
state in the QD 



a\ e''(^+''^"'' Wz 



= -l) 



(5) 



where cr = ±1 are eigenstates of az- 

QD/s-wave superconductor proximity effect The real- 
ization of MPs in the QD system requires the introduc- 
tion of the superconducting proximity effect from neigh- 
boring s-wave SC islands. Let us now consider how the 
proximity effect can by induced from the SC islands in 
the QD/SC chain in Fig.[TJa). Since the states in the SC 
are gapped we can integrate out the states in the SC and 
replace the SC by an effective Hamiltonian , which is 
schematically written as 



i:sc{E^O)^tGsc{E^O)t^ 



(6) 



where i? <C A (A being the superconducting gap in the 
SC island) is the BdG quasiparticle energies we are solv- 
ing for, Gsc is the effective superconductor Green func- 
tion and t is the tunneling between the QDs and the SCs. 
The Green-function of the SC (for uniform A), even in 
the presence of disorder, is written as 



Gscixx ; ^ = 0) = 2^ 



e„ A 
A 



(7) 



where 0„ (x) are the spatial parts of the electronic eigen- 
states with energy e„ relative to the fermi-level in the nor- 
mal state of the SC island and we have used the choice 
{u^{x),ui{x),vi{x),~v^{x))'^ for the Nambu spinor so 
that the BdG Hamiltonian for the s-wave superconduc- 
tor is spin-independent. Here Ua-{x) and Va{x) are the 
electron and hole part of the BdG quasiparticle wave- 
functions. The qualitative form for Gsc is valid even 
for strongly disordered s-wave superconductors and the 
only role of strong disorder in the SC islands is to re- 
duce the coherence length of the superconductor, which 
limits the width of the SC since the SC island is re- 
quired to be shorter than the coherence length so that 
the matrix elements of Ssc oc Gscixx') is non- vanishing 
between neighboring QDs. In particular, for states at 
energies much smaller than the quasiparticle gap A in 
the superconductors, one can assume « 0, so that 
Gsc{E ~ 0) is particle- hole symmetric, Hermitean and 



spin-independent. The effective self-energy induced by 
the superconductor in coordinate (i.e. x,x') space can 
be written as a spin-independent sum of a normal and 
anomalous part as 



^sc{xx') = y!-^^{xx')tx + Y!-g(}{xx')Tz, 



(8) 



where t}-^^^\xx') are the normal and anomalous compo- 
nents of the self-energy. Here the matrices r^^y^z are Pauli 
matrices in the Nambu space containing the particle- 
hole degrees of freedom. To calculate the effective BdG 
Hamiltonian in the sub-space of the low-energy (i.e. near 
the fermi-level ai E ^ 0) active levels in each QD, we note 
that the low-energy basis of the n*'* QD consists of an 
electron-like wave-function described by an x-dependent 
4 component spinor wave-function \E'e_„(a;) = ('!/'n(a;), 0)"^ 
and a hole wave-function "^h.nix) = (0, icTj,';/'* (a;))-^ = 
{0,iayipn{x))'^ (since ipn{x) is real), where ipnix) is the 
position dependent 2-component spinor wave-function of 
the active level in the QD given by Eq. [5] Note that the 
iffy matrix flips the spin of the electron spinor relative to 
the holes. 

Taking the relevant matrix elements of the tunneling 
operator t of the n*^ SC island between the right end 
of the n*^ QD and the left end of the {n + If^ QD we 
find that the effective normal tunneling between QD„ 
and QD„+i is given by 



t>0,n \ ^^SC,n l<^0,ri-|-l) 



7r(a„-ha„+i) 

COS( ; [- Vn ^ Vn+l) (9) 

'so 

where we have used the QD wave-functions from Eq. [5] 
and (00, n I ^scn \4>o,n+i) IS the spatial part of the matrix 
element of the QD wave-function with the superconduct- 
ing self-energy S^^^^. Similarly, the anomalous hopping 
or cross- Andreev reflection amplitude induced between 
the dots is given by 



A„ = (*e,n| Ssc,« |*/i,n+l) ^ S^'^'^ (a„ , -a„+i ) 
. ,7r(a„ -f a„+i) 

Sm( ; \-T]n -??n+l), 

'so 



(10) 



where (0o,n| '^sc n \'t>o.n+i) is the spatial part of the ma- 
trix element of the QD wave-function with the supercon- 

(A) 

ducting self-energy ^scn- limit of weak barrier 

transparency, when the tunneling matrix elements are 
tuned to be of order ^ A, we expect the matrix elements 

(00, n I ^scl'* I0o.ri-i-i) to be of order of the tunneling ma- 
trix elements [13|, which can be made of order A. Note 
that combining the above equation with Eq. [SJ leads to 
the fact that in the limit of vanishing Rashba spin-orbit 
coupling q;, Iso diverges, leading to a vanishing rjn and 
therefore A„. Thus the anomalous hopping amplitude 
A„ is proportional to the spin-orbit coupling. An analo- 
gous calculation leads to the conclusion that the on-site 
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contribution (i.e. from QD„ to QD„) of the anomalous 
self-energy vanishes, while the on-site contribution of the 
normal self-energy just renormalizes /i„ by an amount 
that can be off-set by changing the gate voltage on gate„ 
in Fig. ma). 

Disordered Kitaev chain The resulting effective- 
Hamiltonian for the QD/SC array with a single active 
level in ipnix) in QD„ can be represented conveniently 
as a second-quantized BCS Hamiltonian of the form of 
Eq. [TJ In the effective Hamiltonian H^ff, the sign of 
the hopping parameters i„, A„ associated with SC„, can 
be flipped by a gauge transformation ^''Pm ^o^' 

all m > n. Similarly, i„ and A„ associated with SC„ 
can be interchanged by the particle-hole transformation 
'i/'m — 7> tpm for all m > n, which is accompanied by an 
additional -^,„, t„i ^ -im and -A,„. 

Moreover, even in the case that certain types of spin- 
flip scattering processes (which have not being considered 
here in the previous paragraphs) lead to complex t„, the 
appropriate phase rotations '0m V'me*^" can be used 
to make i„ > 0. Therefore, without loss of generality, 
we can choose a gauge where > |A„| together with 
tn > 0. 

As mentioned earlier, in the perfectly periodic clean 
limit Heff is in the topological phase for < |t„|. 
In this limit, a sign- flip of A^ for all m > n creates a 
TT-junction, which constitutes a domain-wall with a pair 
of localized zero-energy MF modes llj . A configuration 
where A„ has the opposite sign for a block of supercon- 
ducting islands SC^ constitutes a pair of domain-walls, 
which confine a pair of low-energy states whose energy 
goes to zero exponentially in the length of the blocks 
of fiipped signs in A„. Of course, having a gapped su- 
perconducting phase does not guarantee a TS phase. In 
fact, in general determining whether a disordered one- 
dimensional superconductor is gapped requires a calcu- 
lation of the transmission matrix between the ends of 
the one-dimensional wire IJ]. In what follows we will 



consider both analytically and numerically, whether dis- 
ordered Kitaev chains which are free of sign disorder in 
A„ and are otherwise random in terms of A„, tn and /i„ 
support a TS phase with zero-energy MFs localized at 
the ends. We will refer to such Kitaev chains as being 
sign-ordered. 

In the supplementary material, using a transfer ma- 
trix formalism similar to Ref. [14|, we show that the 
Kitaev chain with real parameters considered here has 
a chiral symmetry and therefore can be mapped to the 
problem of the existence of zero modes at the ends of 
a non-Hermitean Hamiltonian with an imaginary vector 
potential 
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In particular, we show that a long sign- 
ordered Kitaev chain has MF zero modes at the end with 
a gap in the bulk whenever sgn(t„A„) = sgn(t„+i A„+i) 
(i.e. sign ordered) and the chemical potential in each QD 
has been tuned using a gate voltage so that /i„ = 0. In 
the case of |/x„| > 0, end modes can be shown to exist 
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FIG. 2: Spectrum of a disordered ID Kitaev chain of length 
= 20 sites. The hopping amplitudes t and pairing am- 
plitudes A are distributed randomly and uniformly in the 
interval [0.5, 1.5] K. The spectrum shows a near-zero-mode 
separated by a gap ~ 1 K in a range of on-site QD chem- 
ical potential < O.lSmey. Thus the TS state shows a 
robust gap in this regime despite the large fluctuations in the 
he 
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FIG. 3: Finite size splitting of Majorana modes as a function 
of chain-length. Finite chain length leads to overlap of end 
MFs which splits them away from zero energy in any finite 
system. The MF splitting Esputung in the disordered ID Ki- 
taev chain (with the same parameters as in Fig. 2) is found 
to go to zero exponentially in chain length A*'. 



whenever 



C = E\og{ 



tn + A, 



-} + Co <o 



(11) 



where the expectation value is taken with respect to 
the random distribution for the magnitudes of |t„|, |A„| 
and Co"^ is the localization length of an effective her- 
mitean Hamiltonian with on-site potential /i„ and hop- 
ping \Jt'^ — A^ . The MF zero mode is localized at the 
ends with a localization length of |C|~^ sites. The param- 
eter regime C < supporting end MFs corresponds to the 
delocalized phase of an effective non-Hermitean Hamilto- 
nian i/o+ A (defined in the supplementary material) [l5{ . 
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The other parameter regime with C > corresponds to 
the locahzed phase 15| of the non-hermitean Hamilto- 



nian Hq + A so that both the non-hermitean Hamiltonian 
Hq and the Kitaev chain H^f f have low-energy modes in 
the bulk. Furthermore, in the supplementary material we 
provide a detailed discussion of why the sign-ordered Ki- 
taev chain with the Hamiltonian in Eq.[l]remains gapped 
for otherwise arbitrary disorder under this condition. 

While the analytic results provide some understanding 
of how zero-energy MF arise in the sign-ordered Kitaev 
chain, it is more insightful to consider the spectra of typ- 
ical realization of finite sign-ordered Kitaev chains. We 
consider in Fig.[2]the spectrum of a 20-site (we have used 
several different number of sites in our numerical work 
getting the same qualitative answer) Kitaev chain with 
open boundary conditions. Choosing a fixed value for the 
gate voltage /i„ across the chain and realistic random val- 
ues of tn and A„, which vary over a significantly larger 
range then their minimum values, we find a robust TS 
phase with end MFs, which appear as near zero-energy 
eigenstates. The end MFs are not strictly localized at 
the ends of the chain but are split by a finite amount due 
to overlap of the decaying MF wave-functions across the 
chain. The exponential dependence on the chain length 
of the energy splitting seen in Fig.[3]shows clearly the ex- 
ponential localization of the MFs and also demonstrates 
that the protection of the topological degeneracy of the 
MFs increases exponentially with the length of the chain. 

The relative signs of the hopping and anomalous tun- 
neling amplitudes i„ and A„ arise the randomness on the 
scale of the fermi wave- vector kp of the electronic wave- 
functions 4>n{x) in the SC islands. Of course, one cannot 
in general expect to eliminate the random relative sign 
arising from such short scale randomness. This prob- 
lem can be remedied by adjusting the individual phases 
of each superconducting island SC„, which according to 
Eqs.|6]and [7]provides an independent handle on the signs 
of tn and A„ . 

Conclusion In conclusion, the two major problems of 
the semiconductor superconductor structures, proposed 
[^, @] in the context of creating non-Abelian solid state 
Majorana modes, are disorder and gating. Both prob- 
lems are solved by using the relatively simple QD/SC 
linear array structure shown in Figure [TJ The use of SC 
grains can potentially introduce random tt junctions in 
the chain, which can produce sub-gap states. These ef- 
fects of strong disorder are avoided by adding an array 
of gates and removing accidental tt junctions by phase 
shifting the superconductors. In addition, the chemical 
potential in each dot is tuned individually with gates, 
thus completely eliminating any gating problem of the 
semiconductor. Finally, structures similar to the QD/SC 
array being proposed in this work are already in use in 
various conventional quantum computing architectures 
such as quantum dot spin qubits and superconducting 
Cooper pair box qubits, thus making this a relatively 



practical (as well as extremely robust) proposal for cre- 
ating non-Abelian modes in the laboratory. 
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Supplementary material 
Zero energy MFs at the end of a ID Kitaev chain 

Consider a ID Kitaev chain Hamiltonian 

(12) 

In this sub-section we prove that if A„ are real and 
sign-ordered i.e. 



previous work 15|. Zero energy (i.e. E = 0) modes 
are obtained as solutions of one of the two decoupled 
equations 



[Ho ± A)C2,i = 0. 



(22) 



sign(A„t„) = sign(A„+it„+i), 



(13) 



and |/x„| < niax(|i„_i|, |A„_i|), then the Kitaev chain 
has zero-energy MFs at the ends. As discussed in the 
main text, we can choose \t„\ > |A„| and i„ > without 
loss of generality. In the topological regime where Eq. [13] 
is satisfied, A„ has the same sign on all sites in the chain. 

Defining a Nambu spinor {ipj^jipn), with a correspond- 
ing particle-hole symmetry operator A = t^K, the BdG 
Hamiltonian corresponding to the above BCS Hamilto- 
nian is written as 



HBdG = HqTz + lATy, 



(14) 



where 



Ho^Yl I'') ^""l + ^"[l'') ('^ + 1| + h.c] (15) 

n 

A = ^A„[|n)(n-M|-/i.c]. (16) 



Applying a unitary transformation U — J^" , HbcIG 
transforms into 



rrrr rj^_( (i?0 + A) 



The BdG equation 

(i/o + A) ^^ ^ Ci \ ^ 



(i?o - A) A ^2 ; ^ V C2 

is now written as 

(7?o + A)C2 = ECi 
(i?o- A)Ci =^C2. 



(17) 



(18) 



(19) 
(20) 



The matrices Hq ± A are written as 
i?o±A = ^-/i„|n) (n| 

n 

+ ^{tn ± A„) |n) (n + 1| + {tn-i T A„_i) \n) {n~l\. 

n 

(21) 

The above Hamiltonians are non-Hermitean and display 
localization transitions similar to the Hamiltonians in 



The sign of A„ can be inverted by flipping the chain, 
therefore in this basis if one of the ends has a MF mode 
with C2 = then the other end has an MF mode with 
Ci = 0. Since both signs of A occur in the above equation, 
we can assume that A„ > without loss of generality. 

Writing the wave-function as "0 = Tim 1*^)' equa- 
tion for the zero-mode can be written as a transfer- 
matrix-like relation 

-^„l/'n + (in+A„)'(/'«+l + (iri-l- A„_i)V'«-l =0. (23) 

Redefining 

Ipn-)- 'ipn= gn'ipn 



9n = 



n 



(24) 

l<m<n 

and dividing Eq. [52]by g„, we obtain the equation as or 
equivalently as 



(25) 

which in turn can be re-written as 

tpn+l 



rj-^ -^W^U - V (iLl " A2_ J^„„l]. 

V - A«) ^ 

(26) 

The above equation represents a zero-mode for a her- 
miteanized tight-binding Hamiltonian 

i?o = ^ -/i„ \n) (n| 



+ ^ ^/{tl-Al) |n) (n + 1| + ^{tl_^ A^i) \n) {n 1| 

(27) 

Let us start by considering the special case of the above 
Hamiltonian Eq. [27l where each gate in the QD has been 
tuned so that fi„ = 0, so that Hq has a chiral symmetry 
which supports zero-end modes similar to Su-schrieffer- 
Heeger model J^]. The equation for ipn then simplifies 
to 



n+l 



(28) 



together with — Q. The above equations can be solved 
for ipp as 



i'2 



(-1)" n /_ _ 



p<n \J (^2p-l 



A^ ^ 



V'2n-1 = 0. 



(29) 
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In the case where the hopping amphtudes are statistically 
independent of each other, the pre-factor 



n 



2p ^2p) 



p<n \/(^2p-l ^ ^2p-l) 



in the large n limit where 



var(log(^(t2-A2))), 



(30) 



(31) 



where the variance var is taken with respect to the ran- 
dom distribution. Since the pre-factor for scales as 
e^, and the factor g„ in the definition of ijjn scales as 
On ~ e~"^^ where 



Ci 



-E\og 



tr, 



tn + Ar: 



(32) 



and E is the expectation value with respect to the dis- 
tribution of tn , A„ , it follows that 



V'2 



-2nCi 



(33) 



so that the zero-mode wave- function V'2n is localized near 
the end of the chain and is normalizable. This shows that, 
as claimed in the main text that localized zero-energy MF 
modes exist for /i„ = 0, for independently ordered bonds. 
However, one should note that such localized MFs don't 
exist for every configuration of t„,A„. For example, a 
dimerized configuration with t2ri-i > i2n would not lead 
to zero modes below a critical value of A„ . 

Now we discuss the more general case for 7^ 0. In 
this case the solution ipn of the transfer matrix relation 
Eq. [21] is expected to scale as 



p"Co 



(34) 



where ^ is the localization length for Eq. [271 Using the 
definition of ^n in terms of V'n we find that 



(35) 



where — {C,i~ Co)j proving Eq. [Tl]so that '0„ is normal- 
ized and localized when ^ < 0. This regime corresponds 
to the delocalized phase of Ref. |l5j • Note that this con- 
dition is also necessary, since in the localized phase one 
will generically obtain low-energy localized states in ffn , 
which will lead to similar energy states in Hq + A [ISj. 
This establishes the conditions where one can rigorously 
expect localized end modes with a gap. However, the 
numerics show that even for short chains and /i„ 7^ 0, 
localized MFs exist. 



Lower bound on the gaps of sign-ordered Kitaev 
chains 



While we have shown the existence of Majorana 
fermions at the ends of the Kitaev chain for a large set 



of parameters, the thermal robustness of the topological 
phase is determined by the quasiparticle gap in the sys- 
tem with periodic boundary conditions. To estimate the 
smallest possible (worst case) gap of Hbcs, we need to 
find a lower bound on to the solutions of Eq. [201 for 
sign-ordered Kitaev chains in the delocalized phase when 
Eq. [TTlis satisfied. 

It follows from Eq. [20l that E^ are eigenvalues of both 
the real-symmetric matrices 

{Ho + Af{Ho + A)C2 - E\2 (36) 

{Ho + A){Ho + AfCi = E\i. (37) 

Since it is easier to prove bounds upper bounds on the 
maximum eigenvalues then lower bounds on the mini- 
mum eigenvalue, we consider the inverse matrices 

{Ho + A)-^^^{Ho + A)-'Ci=E-^Ci (38) 

{Ho + A)-\Ho + A)-^'^C2=E-^C2. (39) 

It is well-known from linear-algebra [l7| . that the maxi- 
mum value of E^^ is related to the maximum row sum 
of the matrix {Ho + A)^^,i.e. 

Cf(ffo + A)-i-^(fro + A)-iCi 



E-^ = 



CfCi 



< max, J2 \iiHo + A)-i),,,|max, ^ \{{Ho + A)-\,\. 

i j 

(40) 

Therefore we get the lower bound on the magnitude of 
the eigenvalues 



\E\ > 



maxj J2 l((^o + A)-i),j|max, Y] \{{Ho + A)-i),^ 







(41) 

Matrices whose inverses decay exponentially and are fi- 
nite have energies that are bounded away from zero. 

The inverse of {Ho + A) can be computed using the 
transfer matrix form for Eq.[251can be written as a matrix 
recursion relation 











1 



1+ ^/(tl-i-^'i-i) 



(42) 

where — {ipm'fpn+i)'^ ■ To calculate the inverse of 
Ho + A for a chain with periodic boundary conditions, 
we can label any site to be and choose 0o = 1 ^^^id 
(j)n^o = 0. Furthermore, we will make the ansatz ^f-i = 
0. This leads to the sequence of equations 







n>0 



(43) 



y/tl-Al 



9 



Using arguments from the theory of ID locahzation using 
transfer matrices 5'„>o is expected to grow as \l/„ ~ e"^", 
where is the locahzation length of the hermitean 
Hamiltonian in Eq. 1271 By transforming back to the orig- 
inal ipn variables using Eq. I24[ we find using Eq. [31] that 
in the delocalized phase 



1 



(44) 



with zeta < so that for long Kitaev chains the correc- 



tions related to our Ansatz =0 are exponentially 
small. Furthermore this proves a bound of ~ on the 
row-sums of max^ \{{Ho + A)~^)ij \ in Eq.HTJ On the 
other hand, by considering the transpose (Hq + A)'^, one 
can obtain a similar bound on the column sum. Substi- 
tuting in Eq. STJ we obtain a gap of order ~ | A„| for the 
Kitaev chain. Of course this is a rough estimate since we 
have not carefully estimated values for the localization 
length Co~^- However, the numerical results in Fig. [Hand 
[31 seem to be consistent with these expectations. 



